arXiv:1503.02609vl [astro-ph.GA] 9 Mar 2015 


Heavy ion irradiation of crystalline water ice 

Cosmic ray amorphisation cross-section and sputtering yield 

E. Dartois 1,2 0. B. Auge 3 ’ 1,2 , P. Boduch 3 , R. Brunetto 12 , M. Chabot 4 ’ 5 , A. Domaracka 
3 , J. J. Ding 3 , O. Kamalou 6 , X. Y.Lv 3 , H. Rothard 3 , E. F. da Silveira 7 , J. C. Thomas 6 


1 CNRS-INSU, Institut d’Astrophysique Spatiale, UMR 8617, 91405 Orsay, France 

2 Universite Paris Sud, Institut d’Astrophysique Spatiale, UMR 8617, batiment 121, 

91405 Orsay, France 

3 Centre de Recherche sur les Ions, les Materiaux et la Photonique (CEA/CNRS/ENSICAEN/Universite 
de Caen-Basse Normandie), UMR 6252, CIMAP - CIRIL - Ganil, Boulevard Henri 

Becquerel, BP 5133, 14070 Caen Cedex 05, France 

4 CNRS-IN2P3, Institut de Physique Nucleaire d’Orsay, UMR8608, 91406 Orsay, 

France 

5 Universite Paris Sud, Institut de Physique Nucleaire d’Orsay, UMR8608, IN2P3- 
CNRS, batiment 103, 91406 Orsay, France 

6 Grand Accelerateur National d’lons Lourds, CEA/DSM-CNRS/IN2P3, Boulevard 
Henri Becquerel, BP 55027, 14076 Caen Cedex 05, France 

7 Departamento de Ffsica, Pontiffcia Universidade Catolica do Rio de Janeiro, Rua 
Marques de Sao Vicente 225, 22451-900, Rio de Janeiro, RJ, Brazil 


Abstract 


Under cosmic irradiation, the ... 
amorphous state. Crystalline ici 
the keV to hundreds of keV ion energies. We experimentally investigate heavy ion ir¬ 
radiation amorphisation of crystalline ice, at high energies closer to tme cosmic rays, 
and explore the water-ice sputtering yield. We irradiated thin crystalline ice films with 
MeV to GeV swift ion beams, produced at the GANIL accelerator. The ice infrared 
spectral evolution as a’ function of fluence is monitored with in-situ infrared spec¬ 
troscopy (induced amorphisation of the initial crystalline state into a compact amor¬ 
phous phase). The crystalline ice amorphisation cross-section is measured in the high 
electronic stopping-power range for different temperatures. At large fluence, the ice 
sputtering is measured on the infrared spectra, and the fitted sputtering-yield depen¬ 
dence, combined with previous measurements, is quadratic over three decades of elec¬ 
tronic stopping power. The final state of cosmic ray irradiation for porous amorphous 
and crystalline ice, as monitored by infrared spectroscopy, is the same, but with a large 
difference in cross-section, hence in time scale in an astrophysical context. The cosmic 
ray water-ice sputtering rates compete with the UV photodesorption yields reported 
in the literature. The prevalence of direct cosmic ray sputtering over cosmic-ray in- 



: interstellar water ice mantles evolve towards a compact 
ice amorphisation was previously monitored mainly in 


* Experiments performed at the Grand Accelerateur National d’lons Lourds (GANIL) Caen, France. Part 
of this work has been financed by the French INSU-CNRS programme “Physique et Chimie du Milieu 
Interstellaire” PCMI) and the ANR IGLIAS 
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duced photons photodesorption may be particularly true for ices strongly bonded to 
the ice mantles surfaces, such as hydrogen-bonded ice structures or more generally 1 
so-called polar ices. 

1 Introduction 



stellar obiects (e.g. Eiroa et al.. 1983t Brooke et all 1996; Dartois et ah. 

998; Brooke et ah. 

1999 


Molinari et ah, 1999t Malfait et ah, 1999; Gibb et ah, 2004; Dartoi 

s, 2005; Boogert et ah. 

2008 


Oberg et ah. 2011, and articles referencing them) and in the warm dust grains 


around evolved late-type OH-IR stars (e.g. 

Soifer et ah, 

198l| Roche & Aitken, 

1984 

Meveretah, 1998; Sylvester et ah. 1999; 

Dartois et ah 

2002; Maldoni et ah. 

2003 

Justtanont et ah, 200a Lombaert et ah, 2013; Suh & Kwoni 20131). 


The interstellar ice mant les present in these environments are imm ersed in a flux of 
cosmic ray particles (e.g. Simpson . 1983 : Webber & Yushakl 1983) that produces new 


1992 

; Moore & FI 

udson, 

2006 

; Fama et ah 

12010 


3) and induces molecules and radicals des- 


2004). 


Inside dense clouds, cosmic rays are the source of the internal FUV (photon energies 
above about 6 eV) photochemistry in addition to radiolysis. Both the FUV photon field 
and cosmic rays colliding with dust grains are major sources of desorption, driving a 
gas to solid state equilibrium, preventing the rapid freeze-out of all the species in the 
solid phase and participating in the construction of molecular complexity. The pho¬ 
todesorption and cosmic ray desorption rates, as well as the corresponding astrophys- 
ical time scales on which these processes are effective, can be derived with dedicated 
ion irradiations of ice films, and the magnitude of the effect can be understood with 
such experimental facilities. 

In this article we investigate the amorphisation of crystalline water ice irradiated by 
swift heavy ions. The experiments performed are described in Sect. 2 and the results 
obtained explained in Sect. 3. These measurements are discussed in Sect. 4 in the 
framework of astrophysical implications. A conclusion is drawn in Sect. 5. 


2 Experiments 

Swift ion irradiation experiments were performed at the heavy-ion accelerator Grand 
Accelerateur National d’lons Lourds (GANIL, Caen, France). Heavy ion projectiles 
were delivered on the IRRSUD, SME, and LISE beam line^J between July 2008 and 
November 2014. These beams were coupled to the CASIMIR (Chambre d’Analyse par 
Spectroscopie Infrarouge des Molecules IRradiees) facility, a high vacuum chamber (a 
few 10 8 mbar) holding an infrared transmitting substrate that can be cryocooled down 
to about 13.6 K. The ice films are produced by placing the cold window substrate in 
front of a deposition line. (Poly-)crystalline ice films were condensed at 150 K on the 


' http://pro.ganil-spiral2.eu/laboratory/experimental-areas 



































































































window, from the vapour phase, then cooled down to the lowest temperature of the 
cold finger and returned to the temperature set for irradiation. We refer to this initial 
state as ’’crystalline ice” within the article. The water ice sample, whose optical depth 
spectra as a function of fluence are shown in Fig. [2] experienced a different deposition 
history. The ice mantle has been formed by annealing of a film deposited at low tem¬ 
perature (13.6 K) of H 2 O co-mixed with N 2 (and traces of CH 4 ) from a previous set 
of experiments. The resulting (poly-)crystalline structure after annealing at 150K - the 
N 2 component sublimates but not H 2 O - and back to 77K is thus slightly different from 
the othe r films, and is ana lysed in the same way. Details of the experimental setup are 


given in Melotl ( 2003 ) and Seneruelo Duarte et al. ( 20091) . The film thickness is chosen 


to give a high band contrast for the infrared absorptions, without saturating the bands. 
(The 40 K experiment infrared signals are close to saturation.) In addition, the chosen 
thicknesses allow the ion beam to pass through the film with an almost constant energy 
loss per unit path length. A Nicolet FTIR spectrometer (Magna 550) with a spectral 
resolution of 1 cm -1 was used to monitor the infrared film transmittance. The chamber 
is equipped with a cryostat that can rotate up to 180 degrees, and the irradiation and 
recording of the transmittance spectra were performed at normal incidence. The pro¬ 
jectiles energies were chosen to cover the high stopping-power range, to complement 
existing data at lower energies with lighter ions. The temperatures were selected to ex¬ 
plore the amorphisation dependencies discussed in the literature for lower projectiles 
energies or electrons. The evolution of the spectra was recorded at several fluences. 



3 Results 


A summary of the experiments performed and considered in the forthcoming discus¬ 
sion is given in Table Q] The column density evolution of the water ice molecules 
submitted to ion irradiation can be described, to first order, as a function of ion fluence 
(F), by a coupled set of differential equations: 


dN/dF = -cr d N - Y s X f 
dvkc/dF = -cr am w c 


( 1 ) 



where N is the total FFO column density, cr d their effective radiolytic destruction cross- 
section (cm 2 ), Y s the H 2 O sputtering yield (FFO/ion), multiplied, to first order, by 
the relative concentration f of water ice molecules with respect to the total number 
of molecules/radicals in the ice film. Here, cr am is the crystalline ice amorphisation 
cross-section,Lind w c the crystalline fraction of the ice. The -<r d N and Y s x f terms 
in the top part of eq. 1 are related to the destruction and sputtering contributions to the 
evolution of the ice column density. Pure water ice is very radiolytic ally resistant. One 
of the main products formed by irradiation is H 2 O 2 , which saturates to a few percent 
of wate r-ic e m olecules, by number, at doses above about 10 eV/molecule (see e.g. 
Moore & Hu dson 200(3 for the production yield of H^O?). The sam e low efficiency is 
observed for FUV irradiation of pure ice bv IGerakines et al. ( 19961) . among others. 

Considering the radiolytic resilience of pure water ice, we impose that f remains 
close to unity. For our considered experimental fluences and ice film thicknesses, if the 
radiolysis is negligible, then the evolution of the water ice spectra is dominated by the 
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crystalline to compact amorphous phase change and sputtering. The 
equations come down to 


to these 


N » N 0 -Y S F 

w c = exp(-cr am F) = N C /(N C + N a ) = N c /N 

where N c and N a are the crystalline and amorphous ice column densities, respectively. 
The column densities of the molecules are followed experimentally in the infrared 
through 



r»v2 

N = ' 

Jv 1 


r(v)dv/A, 


(3) 


where the integral of optical depth (r) is taken over the band frequency range, and 
A is the band strength value (cm/molecule) for the vibrational mode and ice struc¬ 
ture (compact amorphous or crystalline) considered. During the irradiation, the phase 
change from a crystalline to a n amor phous compact state induces a variation in the 
band strengths (see e.g. ' ~~ 

model. 


Leto & Baratta , l2003l) . Therefore, assuming a direct impact 


1 r n 

N = — r(v)dv [w c + cr(l -| w c )]; a = A c /A a 

Ac Jyl 


(4) 


where A c and A am are the band strength values for crystalline and compact amorphous 
ice (integrated over the same wavenumber range). The measured integrated optical 
depth, Iqd = J-j r(v)dv, can then be fitted with 


Iod 


w c + a(l - w c ) 


[No - Y S F]. 


(5) 


To optimise the crystalline to amorphous transition cross section determination, the 
analysis is performed on the spectral region where the contrast between amorphous 
and crystalline ice absorption profiles is maximum. The band strength value variations 
expected over the full profile between crystalline, low-temperature porous amorphous 
and compact am orphous phases are in the several tenths of a percent range (see e.g. 
Leto & Baratta , 2003). However, the most intense profiles changes lie in the range cov¬ 
ering the spec ific crystalline ice peaks, around 3130 and 3210 cm” 1 (see e.g. Fig. 1 in 
Leto & Baratta, 2003, and Fig.[T|in this article). We thus perform the spectral analysis 
over the 3000 to 3250 cm” 1 profile interval (shown in Fig[T}. where the contrast be¬ 
tween crystalline and compact amorphous band strength value variations is maximum. 
Similar contrast approaches are often considered by authors using the near-infrared 
1.65 urn overtone mode (e.g. Leto et all 20051: Mastrapa & Brown , 2006 : Famaetal 
20 HI) . M 


The measured, baseline-corrected, optical depth spectra evolution, in the water ice 
stretching mode spectral region, upon irradiation, are displayed in the upper panels of 
Figs. Q]to[9]for different temperatures and different ions. The corresponding integrated 
optical depth Iod as a function of fluence is shown in the lower panels. The model fit 
is overplotted, and the fit with the amorphisation cross section without the sputtering 
contribution is also shown. The derived cross-sections are given in Table [I] together 
with the sputtering yield estimate when possible. 
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Figure 1: Upper panel: optical depth spectrum evolution upon 46 MeV Ni swift ion 
irradiation at 13.6 K of a thick water ice film crystallised at 150 K. The crystalline ice 
spectrum and the spectrum at the highest fluence are labelled; the intermediate spectra 
correspond to the fluence evolution given in the lower panel. Lower panel: Integrated 
optical depth Iod over the (vi =) 3000 to (V 2 =) 3250 cm~' water-ice stretching-mode 
profile interval, where the amorphous-to-crystalline contrast is higher, as a function of 
fluence (filled circles). A model fit (from Eq. 5) is overplotted (full lines). The fit 
obtained without any water ice sputtering contribution is also shown (dashed lines) for 
comparison. 
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Figure 2: Same as Fig. Q] for 77 K. See the experiments section for the difference in 
deposition history for this sample. 
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Figure 3: Same as Fig. Q] for 100 K. 
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Figure 4: Same as Fig.[l]for 125 K. 
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Figure 5: Same as Fig.[T|for 19.6 MeV Ne swift ion irradiation at 15 K. 
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Figure 6: Same as Fig.[T|for 819 MeV Kr swift ion irradiation at 77 K. 
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Figure 7: Same as Fig.Q]for 3.64 GeV Ni swift ion irradiation at 13.6 K. 
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Figure 8: Same as Fig.Qjfor 3.64 GeV Ni swift ion irradiation at 40 K. 
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Figure 9: Same as Fig.Q]for 81 MeV Ta swift ion irradiation at 13.6 K 
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Table 1: Summary of experiments: E - incident ion energy; Se - electronic stopping power; Sn - nuclear stopping power; T - ice tempera¬ 
ture; cr - compact amorphisation cross-section; K 1 - compact amorphisation dose; Ys - sputtering yield 


Ion 

E 

Se n 

Sn/Se 

IceT 

Ice film 

Projected 

cr 

1/K d 

Y s 

Reference 






thickness* 

range c 




/exp. 


(MeV) 

(eV/A) 

(xl0~ 3 ) 

(K) 

(>m) 

(//m) 

(10 4 A 2 ) 

(eV/H 2 0) 

(H 2 0/ion) 

date 

Ni 26+ 

3640 

62.85 

0.42 

13.6 

0.39 

3.4 xlO 3 

0.16+0.12 

1 25 +3 - 75 

0 54 

_ e 

27-03-2014 





40 

0.62 

3.4 xlO 3 

0.28+0.14 

0 72 +0 - w 

.24 

_ e 

27-03-2014 

Ne b+ 

19.6 

143.1 

1.24 

15 

0.45 

16.1 

0.62+0.09 

0 74 +u - 12 

_ e 

09-06-2012 

Kr ,l + 

819 

409.3 

0.65 

77 

0.22 

175 

1.65+0.26 

0 79 +ii o 

_ e 

11-11-2014 

Ni u+ 

46.0 

461.5 

2.8 

13.6 

0.10 

19.4 

1 .9+1.2 

0.78 + ‘| 3 

2610+1500 

03-05-2014 





77 

0.10 

19.4 

2.65+2.0 

0 56 +172 

-0.24 

0.82^ 32 

1 so 4 ”- 33 

3060+1690 

02-05-2014 





100 

O.U. 1 

19.4 

1.8+0.73 

3200+1080 

03-05-2014 





125 

0.23 

19.4 

1.14+0.23 

9100+2400 

03-05-2014 

Ta 24+ 

81 

793.5 

15.9 

13.6 

0.58 

20.3 

5.02+2.97 

0 51 + ^ 7J 

u - ja -0.19 

17400+11800 

23-09-2011 


a The v alues a re based on SRIM/TRIM calculations for pure ice. The density considered for the crystalline ice is 0.93g/cm 3 
( Feistel & Wagner . 2006)) . * approximate thickness assuming a band strength value of 2.7x 10 1 wn/molecule for the crystalline-ice stretch¬ 
ing mode . c calculated for a hypothetical semi-infinite ice target. d K 1 = Se / (cr n), in eV/molecule e insufficient signal-to-noise or fluence 


to provide a meaningful number 
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Discussion 


4.1 Amorphisation 


The radiation-induced amorphi satio n of crystalline ice has been summarised recently 
for many experiments in Fama et al. ( 2010t) . mainly for light ions (except Ar) and elec¬ 
trons in the 3-800 keV range. The evolution of the compact amorphous phase was 
estimated as a function of irradiation dose D (energy deposited per molecule), follow¬ 
ing an exponential behavior: 


w a 


1 - exp(-K cryst D) = 1 - w c = 1 - exp(-<r am F), 


where 


then 


D[eV/molecule] = 


k: 


F[ions/cm 2 ] x S to tai[eV/cm] 
n[molecule/cm 3 ] 


( 6 ) 


(7) 


( 8 ) 


cryst [ e V/molecule] = S to tal/(c r am x n). 

For simplicity, unless otherwise stated, it is written as K _1 (= K tr [, sl ). In Fig. [lOlwe 
display the amorphisation dose K 1 (in eV/molecule) using our derived cr am , compared 
to p reviousl y repo rted data using the same infrared stretching-mode analysis. 


Fama et al. (12010b point out similarities and discrepancies between results obtained in 


different laboratories, also using various methods to follow the amorphisation process 
(IR spectroscopy, electron diffraction, and channelling). Part of these discrepancies 
are probably attributable to different preparations for the crystalline ice film (annealed, 
deposited at high temperature and lowered to the temperature of irradiation), as well as 
various simplifications for the fitting method to the amorphisation fraction, including a 
maximum fraction of the amorphised ice or not doing so. 

As reported, except for the Strazzulla et al. ( 1992b measurements, these experiments 
yield full amorphisation at high fluences. Our measurements for swift heavy ions are 
in very good agreement with those for light ions at the lowest temperature (13-16 K) 
and with the Ar irradiation at higher temperature. K 1 increases with the temperature, 
showing an increasing resistance to amorphisation with temperature as one approaches 
the amorphous to crystalline transition, but swift heavy ions seem more efficient at 
higher temperatures than do light ions. 

The ion-induced amorphisation of crystalline ice reporte d in the keV to hundreds of 
keV with electrons and light ions (see Fama et al. . 2Q10L summary) shows that an in¬ 
creasing amount of deposited energy is needed to amorphise the ice when the tempera¬ 
ture increases. The heavier Ar ions (100 keV Ar + ; Baragiola et all 2005 ) do not follow 
this trend, which is in line with our results for swift heavy ions, but including a signifi¬ 
cant nuclear stopping-power contribution for Ar. 

Comparing the measured K 1 in the 70-80 K temper ature range, derived for 800 keV 
protons (K ^700, far infrared ice b and s, Moore & Hudson! 1992 ), 100 keV elec¬ 
trons (KT 1 ^ 10, diffrac tion, Lepault et al. . 1983 ). 100 keV Ar + (K 1 =0.7, infrared band, 
20051) . and this work (819 MeV Kr 33+ 


Baragiola et al 


Ni 11+ , K'” 1 =0.56 +1 ' 72 


K- 1 =0.79!° ]5 and 46 MeV 


0 24 ), the cross-section dependence varies by orders of magnitude and 
suggests a more complicated relation than the one related only to the dose, at least for 
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electrons, light ions, and moderate temperatures. K 1 even diverges for electron irra¬ 
diations above 70 K. Such a behaviour also seems true for photons ( Kouchi & Kurodal 
1990), and it has been suggested that the presence of a mechanism that restores the 
original structure explains the strong temperature dependence of electrons and photons. 
One may speculate that defects in the crystalline ice network, induced by the irradia¬ 
tion, may be more efficiently annealed for the lower energy densities deposited by light 
ions and electrons, because the restoration would involve isolated water molecules. 
This would explain the strong temperature dependence of K _1 . This mechanism may 
become much less efficient for high-energy density deposition, such as in the heavy ion 
case, that produces individual tracks of amorphized ice. The investigation of these dif¬ 
ferences should be further explored, because although they are less abundant in cosmic 
rays, if they possess such a flatter temperature dependence, the heavy ions should play a 
more important role at intermediate temperatures for the crystalline ice amorphisation. 
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Figure 10: K amorphisation dose as a function of crystalline ice temperature for the 
present measurements (red filled symbols: triangle, Ni 26+ ; downward triangle, Ne 6+ ; 
cross, Kr 33+ ; circle, Ni 11+ ; square, Ta 24+ ) and previous measurements following the 
evolution of the same ice feature (open symbols: downward triangle, 3 keV He + ; dia¬ 
mond, 30 keV H + ; square,_30 keV He + ; circle, 60 keV Ar 2+ ; upward triangle, 100 keV 
Ar + ; see Table 1 in Fama et akl 20 id) . The purple star shows K 1 for crysta lline ice 


irradiated with FUV Lyman a photons, for comparison (ILeto & BarattaL 2003). 


For temperatures below 100 K, the average crystalline ice K cr L t amorphisatior@ 



* crystalline to amorphous compact ice 
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dose that we deduce is about 0.76+0.24 eV/molecule. The corresponding low-tem peratu re 
porous amorphous ice compaction K,J io| h dose required is 0.24 '{j' 24 eV/molecule (iDartois et al. 
201 3l) . i.e. about three times lower. The corresponding astrophysical time scales to 
evolve from the crystalline to the compact amorphous phase are thus three times larger 
than the on es gi ven for low temperature porous amorphous ice compaction in Table 3 
of Dartois et al. ( 2013 ). 



Figure 11: Water-ice sputtering yield as a function of electronic stopping power. The 
result of our measurements (data above 1000eV/10 l5 H 2 0/cm 2 ) are displayed along 
with the data compiled in Brown et al. (1984, see their Fig. 5 for details on the 
compiled data) at lower stopping power. The fitted dependence of the yield is given 
by the dashed line, and the coloured region shows the standard error of estimate on the 
yield. See text for details. 


In the ion irradiations performed at high electronic stopping power (Ni 11+ and 
Ta 24+ ), the sputtering of ice can be estimated from the infrared spectrum at high flu- 
ence. Th is yield is displayed in Fig. |TT] together with previous measurements reported 
by Brown et al. (19841). 

The data display a power law dependence, which can well be fitted with 



Yield [H 2 0/ion] = aSf [eV/10 15 H 2 O/cm 2 ] 
log 10 (a) = -2.36 ± 0.3;/? = 1.97 ± 0.07. 


(9) 

( 10 ) 
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The sputtering yields for water ice measured in our experiment s are s lightly below 
the yields extrapolated by taking only the data points from Brown et al.l ( 1984b into 
account. This may arise from the differences in ice sample preparation or from the 
fact that the measurements via the infrared are only sensitive to the yield at large flu- 
ence, when the crystalline ice has already been amorphized (into a compact amorphous 
phase) and the first layers have already been sputtered. Part of the difference in the 
yield might thus arise from the ice already having been processed. Such a difference 
could be addressed with direct desorption measurements in the same swift ion range 
with the very first ions impinging the target, using another detection technique. 


4.3 Ion sputtering versus photodesorption 

The sputtering yield dependence can be used to infer the ion-induced sputtering of FEO 
ice in a dense cloud. This is evaluated by integrating the product of this cross section 
with the cosmic ray flux over their energy and abundance distribution. The detail of the 
calculation can be found for the low-temperature porous amorphous ice compaction in 
Dartois et al.l (120131) . so we refer the reader to this article. For the differential cosmic 


ray flux, we use the same functional form. The form parameter Eo(MeV) will affect the 
differential flux of low-energy cosmic rays, with no incidence on the high energy. We 
adopt the galactic cosmic rays (GCR) distribution of cosmic ray abundances. In the 
equations we replace the low-temperature, porous amorphous-ice-compaction cross- 
section by the sputtering yield derived above and derive the FEO ice sputtering rates 
reported in Table [2] 


yield derived above a 

Table 2: H 2 0 -ice spu 


sputtering rates 



J 

GCR 


Eo(MeV) 

200 

400 

600 

^(s-') a 

3.34(-16) 

5.89(-17) 

2.12(-17) 

77 GCR (H 2 0 /cm 2 /s) b 

27.0+ 2 ® 

8 - 3 -4:? 

4 1 +4 1 

^ - 2.1 


a calculated ionisation rate; ' HiO sputtering rate. 


In the shielded part of a dense cloud, the cosmi c rays induced for UV photon flux 
is FUVxTO 4 photons/cm 2 /s (e.g. Shen et al. . 20041) . Therefore, the direct cosmic ray 
water-ice sputtering will contribute equally to a FUV photodesorption yield // FUV for 


^ccr/Fuv 


10 


-3 


(ID 


where ;/ GCR is the CR HiO sputter ing rate. The FUV water-ice photodes orptio n yield 
is still under debate (e.g., Kulikov et al. ( 2010l).<0.05 molecule/ph oton; Oberg et all 
( 20091) . ~ 1, 6 x 1 0 3 molecule/pho ton at T-10 K uWestlev et al] ( 1995 ). ~3-5xl0 3 molecule/photon 
for 35<T<85 K, DeSimone et al. ( 2013 ) ~1 x 10 3 molecule/photon at 100 K for 157 nm 
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photons), because the wavelength-dependent absolute yields and their temperature evo¬ 
lution are more difficult to measure than species such as CO or N 2 , which are less 
bonded to the ice surface. 

In dense clouds, which are shielded from the external interstellar standard radiation 
field (ISRF), the FUV photons are induced by the impact excitation of molec ular hy 


drogen by secondary electrons generated by cosmic rays (e.g. Gredel et al., 198' 


I: 


Therefore, the magnitude of the photodesorption yield and the direct sputtering rate 
by cosmic ray ions are related. From the fitted cross section to the sputtering yield 
measured data, it seems clear that the cosmic-ray-induced water-ice desorption com¬ 
petes with the FUV photodesorption yield. 

5 Conclusion 

Swift heavy ion irradiation of crystalline water ice, at high energies, amorphise the ice 
efficiently and release water molecules by surface sputtering. The main conclusions of 
our study follow. 


- The experimentally measured swift ion dose (K 1 , eV/molecule) required to amor¬ 
phise crystalline ice is about three times higher than the dose for the compaction of 
low-temperature porous amorphous ice, which suggests a better resilience to a phase 
change than the low-temperature porous amorphous ice. 


- The measured temperature dependence of the cross-section for the crystalline to 
compact amorphous-ice transition shows a flatter dependence than measured for elec¬ 
trons and/or light ions. The differences should be investigated further, in particular to 
constrain to what extent the deposited energy density, much higher with heavy ions, 
has an impact on the temperature dependence, in addition to the dose that is often used 
as the reference to compare cross-sections. 

- The water ice sputtering yields, which were derived at large fluence irradiations 
when the ice is amorphized, are in good agreement with measurements at lower elec¬ 
tronic stopping powers using a different analytical technique. The dependence of the 
sputtering yield is quadratic over three orders of magnitude of electronic stopping 
power. 

- Combining, in a simple model, the galactic cosmic ray distribution of ion abun¬ 
dances and the measured water-ice-molecule sputtering yield, the FFO-ice sputtering 
rate in dense clouds is estimated in the range ;/ GCR =4.1 '/ [ H20/cm 2 /s at a relatively 
low ionisation rate *2x 10~ 17 s _1 , and rises to 27.0/ 2 ® H20/cm 2 /s at a higher ionisa¬ 
tion rate ® 3 x KU 16 s~'. 


The magnitude of this sputtering mechanism provides a sputtered FUO-ice molecules 
major feedback on the gas phase, competing with the FUV photodesorption. They are 
nevertheless related, because FUV photons are mainly cosmic-ray-induced in the dense 
regions where water ice is observed. 
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